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Oscillations near an Isosceles-Triangle Solution of the 

Problem of Three Bodies as the Finite 

Masses Become Unequal. 

By Daniel Buchanan. 



§ 1. Introduction. 

In an article entitled " Oscillations near one of the Isosceles-Triangle 
Solutions of the Three-Body Problem," which appeared in the July (1915) num- 
ber of the Proceedings of the London Mathematical Society, the author of the 
present paper discussed the periodic oscillations of an infinitesimal body about 
a straight line drawn through the centre of gravity of two finite bodies of equal 
mass and perpendicular to the plane of their motion, which was assumed to be 
circular. In the problem now under consideration the finite bodies are assumed 
to be of unequal mass and the third body is assumed to be infinitesimal. The 
finite bodies are started so that they move in circles, and the infinitesimal body 
oscillates about the straight line through the centre of mass of the finite bodies 
and perpendicular to the plane of their motion, as in the former article. Initial 
conditions are determined so that the oscillations shall be periodic. The solu- 
tions for the motion of the infinitesimal body are expansible as power series in 
a certain parameter s which represents half the difference in mass between 
the finite bodies. When e=0 the solutions reduce to the simplest case of the 
isosceles-triangle solutions, viz., that in which the finite bodies are of equal 
mass and move in circles and the third body is infinitesimal. 

As we sball have frequent occasion to refer to the former paper we shall 
refer to it as Proc, followed by the number of the equations or the section. 

§ 2. The Differential Equations. 

Let nix an d wi z denote the two finite bodies of mass m — s and m+e, respect- 
ively, and let (i denote the infinitesimal body. Let the unit of mass be so chosen 
that m=l/2, the linear unit so that the distance from m 1 to m 2 sball be unity, 
and the unit of time so that the Gaussian constant shall also be unity. Let the 
6 



42 Buchanan: Oscillations near an Isosceles-Triangle Solution of the 



origin of coordinates be taken at the centre of mass, the plane of motion of 
m x and m 2 as the £»7-plane, and let the coordinates of ni x , m 2 and (i be denoted 
by k\ j Vi > £1 ; £ 2 > % » £2 and I, yj, £, respectively. If the masses m x and m 2 are 
started from the points 1/2+e, 0, 0, and — 1/2 + e, 0, 0, respectively, so that 
they will move in circles, then 



(1) 



£i=(i+e)oos(t—t ), £ 2 = — (i—e) cos (*—£<>), ' 
«7i=(i+«).sia(*--<o)i >7 2 = — (4— e)sin(tf— * ), . 

and the differential equations of motion for the infinitesimal body are 

^-![K^[|+|]KHHH]' 

P» = [(*-*.)'+ O.-^'+S']*. 

If we refer the motion of the system to a set of rectangular coordinates 
rotating about the £-axis with the uniform velocity unity, that is, if we trans- 
form the coordinates of (i by the substitutions 



(2) 



g=as cos (t — t )—y sin (t — 1 ), 
vi=x sin (t — t ) +«/ cos (t — t ), 



(3) 



then, after (1) and (3) are substituted in (2), the differential equations become 

•" = -5R + S +m R-J]' 



These equations admit the Jacobian integral 



(4) 



(x') z + («/') 2 + (3') 2 =a 2 +/-f- — — + iil^- 6 + const. 



(5) 
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When e=0 and the infinitesimal is projected along the g-axis, then x=y==0 
and the differential equations (4) reduce to 

8s 

z" —■ 



(l + 4s 2 )f 

The periodic solution of this equation is, Proc. (4), 

3 

z = ^ — a sin(r — t ) + r^a s [sin3(T' — <t ) — sin(<r— r )J + . . . ., 

(t—to) s 9 , 141 4 , 35 , . 

The constant a is a variable parameter and a/VH(l + ^)i denotes the initial 
projection of (i from the origin when e=0. The numerical values of t and <r 
can both be taken to be zero without loss of generality. As series similar in 
form to 4> occur frequently in the sequel, we shall call them triply odd series 
inasmuch as they contain only odd powers of a and odd functions of odd mul- 
tiples of T. 

§ 3. The Equations of Variation. 

We wish to show the existence of periodic solutions of (4) which are ex- 
pansible as power series in s and which for e = reduce to x = y se 0, z = ^. 
Let us substitute in (4) 



(*-*„) = (T-T„)Vi(l + S), a=*+w, (6) 

where w vanishes with s. If we denote derivatives with respect to t by dots, 
and expand the right members as power series in s, x, y and w, then the 
differential equations (4) become 



x-2Vi(l+S)y-i(l+S)x = x2 Xf>+ 2 2X<'V, 

i=0 i=0 j=l 



y + 2Vi(l + Z)x-i(l + 8)y = y2YP + y 2 2 I?V, 

i=0 **0 /=1 

00 00 oo 

w+W w = 2 W$ 0) + 2 2 W? e j . 

t=2 i=0 }=1 



(7) 



Before characterizing the undefined terms in (7), we shall define a triply 
even series, as such series occur frequently in the sequel. A triply even series 
is a power series in a 2 with sums of cosines of even multiples of t in the 
coefficients. The highest multiple of r in the coefficient of a u is 2 k. Such 
a series is called triply even since it is even in a and only even functions of 
even multiples of t enter the coefficients. 
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The Xf\ Yi 0) and Wf> are homogeneous polynomials of degree i in x, y 
and w, but only even powers of x, and y enter while w may enter to any power. 
The coefficients in Xf> and Yf> are triply even or triply odd series, according 
as i is even or odd, respectively. The coefficients in Wf y are triply odd or 
triply even series, according as i is even or odd, respectively. As far as the 
computations have been carried out we have, Proc. (5), 

Wo = 1— a 2 (|-9 cos 2rj + a*!— -48 cos2r+ ^ cos 4<r) + . . . ., 

Xtf> = 2—3 a 2 (l— 4 cos 2 t) + . . . . , 

7f = -1 _|a 2 (l+2 cos 2*) + . . . ., 

X<» = _48(l + 8) +w, If> = 12(1 + 8) 4 to, 

XP = a 2 [16— 24 a 2 (7— 10 cos 2r) + . . . .]— i/ 2 [24— 36 a 2 (2— 5 cos 2t) + . . . . ] 

— w 2 [24— 108a 2 (4— 5 cos 2v) + .... ], 
r< 0) = — a 2 [24— 36a 2 (2— 5 cos 2t) + . . . . ] +*/ 2 [6— 3 a 2 (1—10 cos 2 *) + .... ] 

+w 2 [6— 9a 2 (7— 10 cos2<r) + ], 

Tff = — a 2 [24a sin<r+ ]+2/ 2 [6a sinr+ ]+w 2 [18a sinr+ ]. 

The Xtf\ rp and WP for j>0 are likewise polynomials of degree i in 
x, y and w, but # enters only to even degrees. If j is odd, the Xtf* are even 
in x, and Yjp and TPp* are odd in x. If j is even, the XP are odd in x, and 
J^ and Wf are even in a;. The coefficients in Xtf* and Y\ n are triply even or 
triply odd series, according as w enters to even or odd degrees, respectively, 
but in Wf> the coefficients are triply odd or triply even series, according as w 
enters to even or odd degrees, respectively. It is observed that X$P = it j 
is even, and Y$> = WP=0 if ; is odd. 

The expansions in the right members in (7) are valid only in a certain 
region of convergence. It is seen from the values of r x and r 2 in equations (4) 
that x, y and z must satisfy the inequalities 

— 1< +s+e 2 -x-2ex+x z +y 2 +z*<l, 
— i<— e+i+x— 2ex+x 2 +y 2 +z*<l. 

The region of convergence is obtained by replacing the inequality signs with 
signs of equality. When referred to the fixed £>?-axes the region is found to 
be the spindle formed by rotating the common portion of two intersecting circles 
about their common chord. These circles have their centres at 1/2 + e, 0, 
and — 1/2+e, 0, and have radii 1/V2. 
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The equations of variation are obtained by putting e=0 in (7) and taking 
only the linear terms in x, y and w in the right members. They are 

£-2Vli(l + S)}y-[*(l + S)+:Xf]*=0, * 

jH-2Vii(i+S)!*-[i(i+a)+rf]j,=o, (8) 

w+W a w=0. 

These equations are the same as Proc. (6) if we put m =i in the latter. 
The solutions of the first two equations are obtained from Proc. (8) by 
putting w =£. They are 

x = A 1 e^-i^ T u 1 +A i e-^ z: i^ T u a +A 8 eP» T u a +A i e~P' T u i , ~i 

where A t (t=l, 2,3,4) are the constants of integration,. u f and v t (*=1, 2,3,4) 
are periodic functions of r having the period 2 n ; and 

&=f/3f*>a 2 * 0=1,2) 

are power series in a 2 with real constant coefficients which are determined by 
the conditions that u t and v t shall be periodic. The values of — ($ 0) ) 2 and 
(/3£ 0) ) 2 are the roots of the quadratic 

64^— 48 y— 119=0, 
and, therefore, 

PF = [ V2-f ]», $°> = [ V2+fJ*. 

The u { and v t are power series in a 2 with sums of sines and cosines of even 
multiples of t in the coefficients, the highest multiple of r in the coefficient of 
a 2 * being 2 k. In u t and v t («=1, 2) the coefficients of the sines are purely 
imaginary and the coefficients of the cosines are real. In u t and v t (i=3,4) 
all the coefficients are real. Further, 

%(V— 1) =m 2 ( — V— 1), u B (r) =Ui(—r), 

«l(V— 1) = «,(- V=T), V,(T) = V 4 (-T), 

«,(0) =l.(t = l,2,3,4). 

The determinant of the solutions (9) and their first derivatives is a con- 
stant,* and at T=0 its value is 



A= 10 V— 119 + terms in a\ (10) 

This determinant is different from zero for a=0 and will remain different 
from zero for a 2 sufficiently small. Hence, the solutions (9) constitute a funda- 
mental set. 

* Moulton's "Periodic Orbits," § 18. 
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The general solution of the last equation of (9) is the same as Proc. (7), 



viz. 



w = A s q> + A,[z + Arty], 
9 

* = COS T + rr; O 2 [COS 3 t — COS T] + 

lb 

g 
X — sinr + — o 2 [sin3r + 5 sinr] + 



(11) 



A 9 2 , I 41 4 , 

where A 5 and A e are the constants of integration. 

§ 4. Existence of Symmetrical Periodic Orbits. 
Let ns choose as initial conditions of (7) 

*(0) = o 1 , A(0) = 0, y(0) = 0, y(0) = a,, w(0) = 0, ti>(0) = «,. (12) 
With these initial conditions it can be shown from the differential equations (7) 
that x is even in r and y and w are odd in <r. Hence, sufficient conditions that 
the solutions shall be periodic with the period 2 n are 

x = y = w = at r = n. (13) 

When the conditions (12) are imposed on the solutions (9) and (11), we find 
that A 1 =A Z , A 8 =A t , A 5 =0, and that A u A s and A 6 are linear functions of 
a;(i=l,2, 3), vanishing with a 4 . Now let us integrate equations (7) as power 
series in e and a t , or, which is more convenient, in e, A x , A s and A 6 . In so far 
as the terms which are linear in Ax , A 3 and A e are concerned, we obtain 

x 10 = ^[eV-i&r Wl +e-V-ift T %] +A s [eP°- T %+e-ft T w«], 
y 10 = V^i^teV-ift'v,— e-V=ifrTt; 8 ]+il 8 [eft'i;,— e-A'v,], " (I 4 ) 

^m = A [#+^ *<*>]• 
When the periodicity conditions (13) are imposed on these solutions, we obtain 

= Ad V^l A+«i(0)] [eV=iA»-e- v-ift»]+j I [/5 f +ii i (0)] [eft*-e-ft»] ' 
+ terms in e and higher degree terms in A x , A s , A 6 and e, 

= V = rv 1 (0)^i[eV-ift»_ e -V-ift»] + v ,(0) J 4 8 [eft*— e-A'] - (15) 

+ higher degree terms in A x , A s , A 6 and e, 

= A 6 An + higher degree terms in A x , A 3 , A 6 and t. 

The determinant of the linear terms in A x , A 3 and A e is 

D = [eV-l/3nr_e-V-l/3i»r] [ e /3 3 T _ e -/S 2 jr] 

{« 8 (0)[V=l/S x +« 1 (0)]-V :=: lt; 1 (0)[/? g +M 8 (0)]}. (16) 
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The last factor has the form 



+ terms in a 2 , 



7V2 

and is different from zero for a 2 sufficiently small. The second factor can not 
vanish as /? 2 is real. Then the determinant D can vanish only when @ x is an 
integer. When /3 X is not an integer, D is not zero; and hence equations (15) 
can he solved uniquely for Aj, , A a and A 6 as power series in s which vanish 
with s and converge for | s | sufficiently small. Therefore symmetrical periodic 
solutions of (7) exist when /? x is not an integer, and they have the form 

00 00 00 

00 = 2 x^, y=l l y j 6 i , w = 2w j e i , (17) 

y=i 3=1 1=1 

where each #,■ , y } and w f is separately periodic with the period 2 n in t. 

From the following considerations of the differential equations (7) we 

show that x and y are odd in e, and w is even in e. The right member of the 

first equation in (7) is odd in x and e, considered together, and even in y. The 

right member of the second equation is even in x and e, considered together, 

and also even in y. Let the solutions (17) be denoted by 

x=x(e), y=zy(s), w=w(e). 

If the signs of x, y and e be changed, the differential equations remain un- 
changed, and therefore, 

x — — x(—e), y= — y( — e), w=w( — s) 

are also solutions. Since the solutions as power series in s are unique, then 

x(e)=— x(— e), y(e)=—y(—e), w{e)-w{—e), 

or, x and y are odd in e and w is even in e. Hence, the periodic solutions of 
(7) have the form 

00 00 00 

x = 2 x 2j+1 e 2 >'+\ y=Xy 2j+1 e^ +1 , w = Xw 2j ^. (18) 

j=o j=o i-x 

The same result would be obtained in the construction of the solutions if the 

forms (17) were used instead of (18). 

§ 5. Existence of Symmetrical Periodic Orbits when /3 a is an Integer. 

When /?i is an integer, the determinant (16) vanishes and the 'terms in 
(15) of higher degree in s, A x , A s and A 6 must be considered in order to 
establish the existence of symmetrical periodic orbits. We take the same 
initial conditions (12), and as in the previous section we integrate (7) as power 
series in A lf A z , A 6 and e. In addition to the linear terms already obtained 
in (14), we need the quadratic terms in A 1} A 3 and A 6 , and the term in A\ 
in so far as it enters x. If these terms are denoted by x u , y u , w n and a% 
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respectively, then all these terms except w n can be obtained from Proc. § 6 by 
putting w =l/2, 1*=0 in the solutions there denoted by the same notation. 
Thus, 

X n = -r J. 1 ^ 6 Q 1 T[eV-lftr^ l _ e _V-lftT M2 ] 

+ — x— A 3 A 9 Q 2 t[ eP* T u 3 —e-^ T u i ], 



y 11 = Xl AlA6 Q lT[e V-lfrT Vl + e -^-if tl r V2] 



+ - 1 r- A z A 6 Q 2 t[e^ v s + e-^ v,], 



(19) 



x 30 = i^?&*[e v - lftT % _ e -V-iftr Uf]t 

where Q lt Q 2 and Q 3 are power series in a 2 with real constant coefficients. 
The term w n is the solution of the differential equation 

Wn+W w u =W n , (20) 

where the right member is a linear function of x\ , yl , wf multiplied by a 
triply odd series. The complementary function of (20) is the same as the 
solution of the last equation of (8), and on using the method of the variation 
of parameters we obtain 

Wi^rlAltyi+Altygj+Alinon-iperio&ic terms j + periodic terms, 

where $ x an ^ $3 are P ow er series similar to $ in (11). 

When the periodicity conditions (13) are imposed on the solutions 

# = 210+351! + %)+ , 

2/= */io +2/U+ , 

we obtain the equations 

0=^,[ft+«,(0)][eft»-e-A*+^^ii,g J [eft*+e-ft*]j 

i-^AtV^lfr+M^O)] [A 6 Q 1 +A\Q 3 ] + terms in s 
and cubic and higher degree terms A-y , A 3 and A & , 



= V— lv 3 (0)A 3 L^— e-ft*+£-4 8 &[eft* + e-0i*]j 

+ ^-^A[±2v 1 (0) A 1 Q 1 +v 3 (0) A,Q a [e**+e-**] 

+ terms in e and cubic and higher degree terms in A x , 
A s and A 6 , 
= n[A A 6 + Alfain) +Al$ 3 (n) ]+ terms in e and higher degree 
terms in A x , A 3 and A 6 . 



(21) 
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Where the double sign occurs the + is to be taken if & is an even integer, 
and the — if it is an odd integer. 

The last equation of (21) can be solved for A 6 as a power series in A\, 
A\ and s which vanishes with A x and A 6 but not with e. "We shall refer to this 
solution as (21c). When it is substituted in the second equation of (21) we 
obtain an equation in A 1} A s and e which vanishes with Ax and A s but not with 
s, the terms of lowest degree being A\ and A s . As the coefficient of A z in this 
equation is a power series in a 2 with additional terms in 1/a 2 , it will be different 
from zero, in general, and the equation can be solved for A 3 as a power series 
in A\ and e which vanishes with A x but not with s. Denote this series for A 3 
by (21b). After substituting (21c) and (21b) in the first equation of (21), 
we obtain an equation (21a) in Ax and s which vanishes with these terms and 
in which the lowest power of Ax alone is A\ • The coefficient of A\ in this 
equation is a power series in o 2 with additional terms in 1/a 2 and will, in 
general, be different from zero. Hence, (21a) can be solved for A x as a power 
series in f I which vanishes with e and converges for | e | sufficiently small. 
There are three solutions for A lf but only one is real, the other two being 
complex. When this power series for A x is substituted in (21b) and (21c), 
A s and A 6 are likewise power series in si . Therefore, when & is an integer ? 
periodic solutions exist having the form 

x = I x li) e*, y=X y^ <? T , w = I w (i) s* ( 22 ) 

/=i /=i j=i 

where each x (p , y ij) and w iiy is separately periodic with the period 2n in t. 

These solutions converge for s sufficiently small numerically. There is only 

one set of real solutions, the other two sets being complex. 

In the practical construction of the solutions it can be shown that x and y 

are odd in si, and that w is even in si . This property of the solutions can be 

deduced directly from the differential equations (7) by an argument similar 

to that used at the end of § 4 if we replace s in that section by si . Then we 

may consider the solutions to have the form 

x = 2 x (2i+1) r 5- , y = 2 y (2i+v s~^~, w = 2 w (2j) s s . (23) 

§ 6. Proof that all the Periodic Orbits are Symmetrical. 

We shall now consider the existence of periodic orbits in which the infini- 
tesimal body is projected from the xy-plane in a direction which may not be 
perpendicular to. the tc-axis, and from a point which may not lie on the #-axis. 
These orbits, if they exist, are called the general orbits. They have the initial 
values 

«(0) = O!, x(0) = a 2 , y(0) = o 9 , y(0) = a iy w(0) = 0, w(0) = a s . (24) 



(25) 



^ ] (27) 
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The initial value of w can be taken to be zero without loss of geometric 
generality since the infinitesimal body must cross the a;t/-plane if the motion 
is to be periodic, and hence the initial time can be chosen as the time when the 
infinitesimal body crosses the #«/-plane. Sufficient conditions that solutions 
having the initial values (24) shall be periodic are 

x(2n) — #(0)=0, y(2n) — */(0)=0, w(2n)—tv(0)=0, 
x(2n)—x(0)=0, y(2n)—y(0)=0, w(2n)—w(0)=0. 

We shall now show that one of these conditions, viz., w{2n) — w(0) =0, can be 
suppressed. 

So far no use has been made of the integral (5). It is by means of this 
integral that we show that one of the conditions in (25) is redundant. When 
(5) is transformed by the substitutions (6), it takes the form 

a?+y 2 +(i>+w) 2 +F(x,y 2 ,w; e)=0, (26) 

where F is a power series in x, y i , w and e having periodic coefficients. Let 
us make in (26) the usual substitutions 

x=x(Q)+x, y—y(0)-\-.y, w=0+w, 

x-x(0)+^, y=y(0)+ij, w=w(0)+l 
where x, x, y, y, w and w vanish at t = 0. We shall denote the equation 
resulting from substituting (27) in (26) by (26a). By putting t = we 
obtain from (26a) an equation (26b) connecting the terms of (26a) which 
are independent of x, ...., w. On substituting (26b) in (26a) we obtain 
an equation 

G(x, t, %, y~, w, w) =0, (28) 

in which at v=0 or 2n there are no terms independent of the arguments 
indicated. The coefficient of w(2n) in this equation is 2[a-\-iv(0)], and it 
vanishes only when *i;(0) = — a. If w (0)-— — a then i(0) =0, since z=^-\-w 
and 4 , (0)=a, and as s(0)=0 it follows that s^0. As we desire solutions 
which are not identically zero we must take w (0) = — a. Therefore the coefficient 
of w{2n) in (28) is not zero at r=2n, and by the theory of implicit functions 
this equation can be solved uniquely for w{2n) as a power series 

W(2n)=H\x(2n), x~(2n), y(2n), y~(2n), w(2n)\, (29) 

which vanishes with the arguments indicated and converges for sufficiently 
small values of the moduli of these arguments. When the first five conditions 
of (25) are satisfied it follows from (29) that the arguments of H are all zero, 
and therefore w(2n) = 0. Then, since w(2n) =iv( 2 n) — w(0), it follows that 
the condition w(2n) — m>(0)=0 is a consequence of the other conditions in 
(25) and may be suppressed. 
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We shall now integrate equations (7) as power series in a t (*=1, . . . ., 5) 
and s. The differential equations from which the linear terms are obtained 
are the same as the equations of variation, and consequently the solutions are 
the same as (9) and (11). When the initial conditions (24) are imposed on 
these solutions, we find that A & =0 and that the remaining A { are linear func- 
tions of <Xj which vanish with a 4 . Then, as in § 4, it is more convenient to 
integrate (7) as power series in A t and e. On imposing the necessary condi- 
tions of (25) upon (9) and (11) we obtain 

= ^,[e«V=iA*— 1] +^ 2 [e-2V-i/3 1 7T_i] 4.^ 3 [ e 2ft7r_i] 

+A i le-***—l'\ + ...., 

= [V :i l/3 1 +M 1 (0)]{^ 1 [e2V- 1 AT_l]— 4 2 [e-2V-i&T— l]} 

+ W2V 8 (0)+v a {0)]{A 3 [e*P.«-l]-A l [e-*P**-l\}+ . . . ., 
= V =: iv 1 (0){^ 1 [e2V- 1 /3 lff _i]_^ 2 [ e - 2 v-iAT— 1]} j. (30) 

+ v 9 {0){At[&**—l]—A t le-*k'-l]}+.. . ., 

= [— /9 1 v 1 (0) + V^tf 1 (0)]{4 1 [e«V-ift*_l]+^[e-«V=ift*--l] j 

+ W 2 M0)+v s (0)]{A s [e^--l]+A l [e-^-l]}+...., 

Q = 2nAA,i>(0) + 

The determinant of the coefficients of A t is 

2n&Ai>(0) [ e 2Vri/3 lff _i] [ e -«V=ift*_i] [e«ft»_i] [ e -«A»_ 1], (31) 

and it vanishes only when /3 X is an integer since /3 2 is real and A, A and $(0) are 
different from zero for a not zero but sufficiently small numerically. If & is 
not an integer the determinant (31) is different from zero and equations (30) 
can be solved uniquely for A i as power series in e which vanish with e and 
converge for |e| sufficiently small. Hence, when & is not an integer, the 
general orbits exist uniquely and have the same form as the symmetrical orbits. 
Since the general orbits include the symmetrical and both classes of orbits are 
unique, then all the periodic orbits are symmetrical when & is not an integer. 
When & is an integer the coefficients of A x and A z in the first four 
equations of (30) vanish and in order to discuss the existence of general 
orbits, it is necessary to consider the terms of higher degree in (30). In 
addition to the linear terms (9) we require, as in § 5, the quadratic terms in 
Ai(i=l, 2, 3, 4, 6) which contain t as a factor, and the terms of x containing r 
as a factor and of the third degree in A± and A 2 . These terms are found in 
the same way as the corresponding terms were determined in §5. They have 
the form 
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% = ^ i 6 Ci ^ [^i e vri ^ T «i-i 2 e-V-i ft r Mj] 

.— — ,. 

+ — ^—A 6 Q % <t [A, eh T u,—A t e-fc T w 4 ] , 



y 11 = ^—±A 6 Q 1 r[A 1 e^^v 1 +A 2 e-^-i^v 2 ] 



(32) 



% = r[Ax AzQi+AsA^s] +A% [non-periodic terms] 

+ periodic terms, 

Xso = -^A 1 A 2 Q s t[A 1 e^-^^u 1 — ^eV-ift^]. 

When the necessary periodicity conditions of (25) are imposed on equations 
(32), (9) and (11), we obtain 

= A a [0**-1] +A[e-tft*-l] + ^i,Ci[i,-i,] 

-f terms in e and higher degree terms in A if 
0=[ft+«,(0)]{^,[e«ft»-l]-il 4 [e-iA*-l]} 

+ ^l [V^^+m^o) j [i.ft+^i.e,] Ui+^,] 

+ ?^El?^g t [ / 3,+« t (0)][ J 4 B e«Ai*+4 4 e-iP.-] 

+ terms in e and higher degree terms in A { , 
= Vi(0)^ g [ 6 «A*--l]— ^Je-«A»— 1]} 

+ 2VEI*^ ( oM i g l u 1 +^] ^ (33) 

+ ^=^v,(0)^,g jl U,e«ffi*+^4e- 1 * ,r ] 
+ terms in e and higher degree terms in A t , 
0=[ftt;,(0)+t;,(0)]{4,[e«ft»-l]+4 4 [e-«A'-l]} 

+ ?^^0 1 [V=iM(0)+%(0)][A-A] 

+ terms in e and higher degree terms in A t , 
= 2nAA s i>(0)+27t[A 1 A 2 (}> 1 (0)+A s A i <l> 3 (0)] 
+ terms in e and higher degree terms in A t . 



Problem of Three Bodies as the Finite Masses Become Unequal. 53 

The last equation of (33) can be solved for A 6 as a power series in A x , A 2 , A 3 , 
Ai and s, and the terms of lowest degree in A 4 are A x A z and A 8 A± . The deter- 
minant of the coefficients of A s and A± in the third and fourth equations is 

2»,(0)[ftv 8 (0)+t>,(0)]Le«A'-l][e-«ft'-l], 

and it is different from zero for a 2 sufficiently small since /3 2 is real. Hence 
these two equations, when A 6 has been eliminated by the last equation, can be 
solved for A z and A 4 as power series in A x , A 2 and s in which the terms of lowest 
degree in A x and A 2 are A\A % and A X A\. When these solutions for A 3 , A t 
and A 6 are substituted in the second equation, the terms independent of e con- 
tain A X A 2 or A X A\ as a factor. Hence this equation can be solved for A 2 as 
a power series in A x and s, and the term of lowest degree in A x is linear in A x . 
Finally, when the solutions for A z , A s , A t A e are substituted in the first equa- 
tion, we obtain an equation in A x and s in which the terms of lowest degree in 
A x is A\. As in (25a) the coefficient of A\ in this equation is a power series 
in a 2 with additional terms in 1/a 2 and will, in general, be different from zero ; 
therefore this equation can be solved for A\ as a power series in ei which 
vanishes with s and converges for |e| sufficiently small. There are three 
solutions as in the symmetrical orbits, but only one solution is real, the other 
two being complex. Since the real solutions for the general orbits and the 
symmetrical orbits are both unique and since the general include the sym- 
metrical orbits, all the periodic orbits are symmetrical when & is an integer. 
Consequently all the periodic orbits are symmetrical whether /3 a is an integer 
or not. 

§ 7. Constructions of the Solutions when /3 X is not an Integer. 

Let us substitute (18) in (7) and equate the coefficients of the various 
powers of s. We obtain a series of differential equations which can be inte- 
grated step by step and the constants of integration arising at each step can 
be determined, as we shall show, so that the solutions shall be periodic and 
shall satisfy the symmetrical initial conditions 

%=y=:w=0 at r=0. 
When these initial conditions are imposed on the solutions (18) we obtain 
^ + i(0)=!fc + i(0)=0, 0=0, ....oo), 1 

w„(0)=0, = 1 oo). J {64c) 

The differential equations for the terms in e are 



«i-2V i(l+a)y i-[i(l+8)+Xf]g l r=Xf> > 
yi+2Vi(l+S)x x -[i(l+8)+Ynyi=0, 
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where Xo° is the same even power series as that represented in (7) by the 
same notation. The complementary functions of (35) are the same as (9), and 
the particular integrals can be obtained by the method of the variation of 
parameters.. The complete solutions are thus found to be 

x 1 =Ai i:> e^~^^u 1 +Ai iy e-^~^n 2 +Ai 1) e^u 3 +Ai 1) e-^u i +G 1 (r) , 1 

where -4|°(i=l, 2, 3, 4) are the constants of integration, and C x (t) and 8 x (t) 
are triply even and triply odd series, respectively. As such series occur 
frequently in the construction, we shall denote them by Cj(r) and 8 f (<t), 
respectively. Since & is not an integer in this case and 8 2 is real, the terms 
of the complementary function must be suppressed in (35) by choosing AP=0, 
(»=1, 2, 3, 4), and therefore the desired solutions become 

* 1 =C 1 (t), ^=8^). (36) 

The differential equation for the terms in e a is 

w 2 +W l) w 2 =W 2 , 

where W t is a triply odd series. The complete solution of this equation is 

obtained by employing the method of the variation of parameters, and it is 

found to be 

w 2 =£f>$+£f> [x+A^-m^+S^r), (37) 

where Bf* and 2?f } are constants of integration, and p 2 is a power series in a 2 
with real constant coefficients. In order to satisfy the periodicity and the 
initial conditions (34), the constants of integration must have the values 

J5f>=0, JSf>=^, (3,8) 

in which case the solution (37) becomes 

W2=4&(r), (39) 

where 8 t (v) is a triply odd power series. 

The remaining steps of the integration can be carried on in the same way. 
The solutions for x n _ x and y 2 j-i are obtained from the particular integrals 
alone since /? 2 is real and 8 t is assumed to be a number which is not an integer. 
These particular integrals are triply even and triply odd series, respectively, 
except for a factor ar 2i >~ 2) which is introduced through the factor a~ 2 in (39). 
The general solution for the w n is similar to (37) and it can be made periodic 
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by a proper choice of the constants of integration as in (38). At the steps 
2j — 1 and 2/ the solutions are 

1 



X *l-1 — ^2j=4 ^W-l ( T ) » 

2/2,-1— ' a u=T ^i-iir) , 



(40) 



§ 8. Construction of the Solutions when /? x is an Integer. 

Let us substitute equations (23) in (7) and equate the coefficients of the 
various powers of ei As in the previous section we obtain a series of differ- 
ential equations which can be integrated step by step, and the constants of 
integration can be determined, as we shall show, so that the solutions shall be 
periodic and shall satisfy the symmetrical initial conditions 

x(Q)=y(0)=w(0)=0. 

When these conditions are imposed upon (23) we obtain 

i«'+»(0)=y«'+»(0)=0, 0=0, ....oo), 1 

w^(0)=0, (j=l, ....oo). J { } 

The differential equations for the terms in ei are the same as the first 
two equations of (8) if we use the superscript 1 on » and y. Since & is 
assumed to be an integer, the solutions of these equations which are periodic 
and which satisfy (41) are 

t/<i> = V-lA (1) [e v - Jft X-e~ V ~ lftT t> 2 ], J 

where A^ is a constant which is undetermined at this step. The solutions 
(42) are real if A{^ is real. 

The differential equation for the terms in el is 

w m + W w a > = W ai =(Ai l >) , S m (r). (43) 

As functions similar to S (2) (r) occur frequently in the construction of the 
remaining solutions, we shall define the function S 2i (r) to be a homogeneous 
polynomial of degree 2i in the exponentials e^ - ^' 7 and e~^~ ie,T of which the 
coefficients are power series in odd powers of a with sums of sines and V — 1 
times cosines of odd multiples of r in the coefficients. The highest multiple of 
r in the coefficient of « 2 * +1 is 2k+l. The coefficients of [e v ~ ]ftT ]'' [e-^-'A*]* and 
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[ e V-i&Tj*|- e -V-iftT-] }} ^+ ^ =2*, differ only in the sign of V^. The part of 
S m (r) which is independent of the exponentials is a triply odd power series. 
If the exponentials are expressed in trigonometric form, then 

S <2i) (<r) = 2 2 2 a^Sffi, t sin i 2fcft ± (21 + 1 ) } *, 

/=0 1=0 ft=0 

where StfX t are real constants : that is, S (2iy (t) is the same as a triply odd 
power series except that the highest multiple of r in the coefficient of a 2/+1 is 
not 2i+l but 2(i+ift) +1. 

The general solution of (43) is found in precisely the same way as (37) 
was obtained. It has the form 

w^=B^+B^[x+Ar^]-a(A^)Y 2) ^+(A^rS^(t), (44) 

where p i2) is a power series in a 2 with real constant coefficients and Si 2) (r) is 
similar to <S Y(2) (t) . On imposing the periodicity and the initial conditions (41) , 
we have 

B?> = 0, Bf=^j{A[») 2 p i2) . (45) 

Then the desired solution of (43) takes the form 

w m = -^(APyS' 2) (r), 

where 8 (2) (t) is similar to 8 (2) (<r) . Thus Ai l) remains undetermined at this 
step also, and the terms in e 3 must be considered. 

The left members of the differential equations for the terms in f 3 are the 
same as the first two equations in (8) if we use the superscript 3 on a; and y. 
Let us denote the right members by X (3) and Y (3 \ respectively. Then 

Z (3) =(^ 1 (1) ) s p (3) +X< 1) , r (3) = V :=: T(^ 1 (1) ) 3 (T <8) . (46) 

The function X$ is the same triply even power series as that denoted in 
(7) by the same notation. The functions p (2i+n and <j <2 ' +1) are homogeneous 
polynomials of degree 2i+l in the exponentials e y/ ~ 10lT and e -V ~ t/3iT , and the 
coefficients are power series similar in form to m x and u 2 in (9). In p (2 ' +1) the 
coefficients of [ e vr "^ T ]>l>~ V ~ lftT ]* and [e^'*]*^" 1 * 1 ']', ]+h=2i+l, differ 
only in the sign of V^ ; while in a (2i+1) the coefficients of [ e ' / - l ^[—e-" / -^ T ] h 
and [e v ~ lAT ]*[— e- v=i * T ] / differ only in the sign of V^l. 

The complementary functions of the differential equations in x iS) and t/ (3) 
are the same as (9). Since X (8) and Y m contain terms which have exactly the 
same period as the periodic parts of the complementary functions, the par- 



(47) 
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ticular integrals will contain non-periodic terms. By the method of the varia- 
tion of parameters we find that the general solutions for x m and y (S) are 

x m =A^e y/ ~ l ^ T u 1 +Af >e- V:ri * T M. 2 + A? e^u s +Af >e-^ T M 4 

+t[ (4P) 3 P?>+Pf J [e v - 1 ^ T w 1 -e- v -^ T w 2 3 +p (S) , 
^<s) = V^I [A?> e v - 101 X-^f e- v ~ 1AT v 2 ] +^4f >e^v 3 - Jf >e-^ 4 
+ V^-rt (^F^Pf+PfJ [e vri01 X+ e- v - 1Ar v 2 ] + V^ff' 3 ', 

where ^| 8> (t=l, 2, 3, 4) are the constants of integration, Pf> (*=1, 2) are 
power series in a 2 with constant coefficients, and p <3> and 5 (3) are similar to p <3) 
and <r (3> , respectively. In order that the solutions shall be periodic, then 

Ap=AP =0, P?> (4«) 3 +Pf =0. 
The solutions for A[ iy are 

AP=P?\ oPf or o 2 PP, 

where Pf > is a power series in a 2 with real constant coefficients, and o, g>* are 
the imaginary cube roots of unity. As the imaginary solutions for A^ lead 
to imaginary orbits, we retain only the real solution. Then, on imposing the 
initial conditions (41) on (47) we find that A¥ l) =AJP, and the desired solutions 
for a? (3) and y (8) become 

The constant Af> remains arbitrary at this step and the terms up to ei must 
be considered before it can be determined. The solutions for # (3) and y &) are 
real provided Af* is real, which is found to be the case. 

The general solutions of the differential equation in el is obtained in the 
same way as (44) was found, and it can be made to satisfy the periodicity and 
initial conditions by a choice of the constants of integration as in (45). When 
these conditions have been imposed, the solution for w u) takes the form 

a 

The differential equations in el are the same as (8) in the left members 
if we use the superscripts 5. If we denote the right members by X (5> and Y a \ 
respectively, then 

where p£ 3) , p£ 6) and oj?\ a£ 8) are similar to p (2 * +1) and <r <2fc+1) , respectively. The 
8 
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complete solutions of the equations in x i5) and y i6) are readily found in the 
same way as (47) were obtained. They are: 

x m =A?h^-^ T u 1 +A^e- v -^ T u 2 +Ate^ T u a +Afe-^ T u i 



ui a ' 



^^ V_l [A^e^-^v.-A^e-^-^v,] -Mfe feT v 3 -,4f e-^< 



V^T* [ p(5)^( S ) + P (5) ] [e V-i A r v +e -V-i A r ] + V; 

ar e 



(48) 



where P| 6) (t = l, 2) are power series in a 2 with real constant coefficients. The 
solutions for x ib) and t/ (6) can be made periodic by putting A^ ) =Af ) =0, and 
then choosing Af> so that p(5)^(3) , p(6)_o 

from which it follows that Ai 3) is real. From the initial conditions (41 ) we 
have A^=A$K 

The remaining steps of the integration can be carried on in the same way. 
The solution at the step 2) can be made to satisfy the periodicity and the 
initial conditions by a proper choice of the constants of integration arising at 
that step. The solution for w 2j contains the constant Af t ~ 1 ' 1 which is not deter- 
mined until the step 2j-\-l. The non-periodic parts of the complete solutions 
at the step 2j-\-l are similar to those in (48), and A ( ? i ~ v enters these solutions 
with the same coefficient as Af> enters (48) except for the factor l/a 2, '~ 2 . The 
solutions can be made to satisfy the initial and periodicity conditions by 
putting ^(2y+D_^(2j+i) A &i+l '*=A 9S+u =Q 

and then determine the Af'~ l) so that the coefficient of r shall be zero. The 
desired solutions at the steps 2j and 2j-\-l are 



a .(2, + l) = ^(2,+l) [e V-lftr Mi + e -V-lAr M2] + ___ 

a 



(2i+n 



^-l&'-ji. — p-V-lftT,,, I I Y l ff (2;+l) 



(2,- + l) = y_ 1 ^p+l)|- e V- 1 Ar g-V-lftr -] + V— -L g 

a 2j-2 

Thus the integration can be carried on to any desired degree of accuracy. 

The solutions of (7), whether as power series in e or si, must satisfy the 
integral (5) identically in s or si. Thus this integral, besides being of use in 
proving that all the periodic orbits are symmetrical, serves as a check on the 
computations. 

Queen's University, Kingston, Ontario, November 15, 1915, 



